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Sub-surface Oxygen and Surface Oxide Formation at Ag(lll): 
A Density-functional Theory Investigation 
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To help provide insight into the remarkable catalytic behavior ol the oxygen/silver system for het- 
erogeneous oxidation reactions, purely sub-surface oxygen, and structures involving both on-surface 
and sub-surface oxygen, as well as oxide-like structures at the Ag(lll) surface have been studied 
for a wide range of coverages and adsorption sites using density-functional theory. Adsorption on 
the surface in fee sites is energetically favorable for low coverages, while for higher coverage a thin 
surface-oxide structure is energetically favorable. This structure has been proposed to correspond 
to the experimentally observed (4 x 4) phase. With increasing O concentrations, thicker oxide-like 
structures resembling compressed Ag20(lll) surfaces are energetically favored. Due to the relatively 
low thermal stability of these structures, and the very low sticking probability of O2 at Ag(lll), 
their formation and observation may require the use of atomic oxygen (or ozone, O3) and low tem- 
peratures. We also investigate diffusion of O into the sub-surface region at low coverage (0.11 ML), 
and the effect of surface Ag vacancies in the adsorption of atomic oxygen and ozone-like species. 
The present studies, together with our earlier investigations of on-surface and surface-substitutional 
adsorption, provide a comprehensive picture of the behavior and chemical nature of the interaction 
of oxygen and Ag(lll), as well as of the initial stages of oxide formation. Copyright (2003) by The 
American Physical Society. 



PACS Nrs.: 68.43.-h; 82.65.-|-r 
I. INTRODUCTION 

Investigation of the interaction of oxygen and silver, 
with regard to the role it plays in technologically impor- 
tant heterogeneous catalytic reactions such as ethylene 
epoxidation and the partial oxidization of methanol to 
formaldehyde, has a long history.^ Despite considerable 
efforts, the 0/Ag system is still not well understood on 
a microscopic level, and the identification of active oxy- 
gen species involved in the above-mentioned reactions, 
remains unclear. Of particular interest is to understand 
the role and nature of the elusive sub-surface O species 
in the function of silver as an oxidation catalyst. Sub- 
surface oxygen is thought to play a crucial role in the 
epoxidation of ethylene (i.e., formation of ethylene oxide, 
alternatively known as epoxide) , which takes place under 
atmospheric pressure and temperatures of 500-600 K, and 
in the formation of formaldehyde from methanol, which 
can occur either by direct dehydrogenation or by oxi- 
dehydrogenation, and also takes place under atmospheric 
pressure and temperatures of 700-900 K and 800-900 K, 
respectively. 

Recent studies of oxidation reactions over certain other 
metal surfaces have found that metal oxides form and 
actuate the reaction, in contrast to the hitherto believed 
pure metal, i.e., for the oxidation of carbon monoxide 
over Ru.^^'"'^^ One may therefore wonder if similarly sil- 
ver oxides form and play an important role. The most 
stable oxide of silver, Ag20, however, reportedly ther- 
mally decomposes at 460 K and atmospheric pressure, 
making this seem unlikely; ruthenium dioxide on the 
other hand is stable up to temperatures greater than 



1000 K.^^ To date there have been a limited number of 
ah initio studies of sub-surface oxygen adsorption and 
oxide formation at transition metal surfaces. These in- 
clude the investigation of oxygen incorporation into the 
Ru(OOOl) surface^^'19 which predicted that after 1 ML of 
oxygen on the surface, O adsorbed in sub-surface sites 
below the first Ru layer where an attractive interaction 
between the oxygen atoms was identified, suggesting is- 
land formation. Within these structures, the overall O- 
Ru bondstrength was reported to be significantly weaker 
than pure on-surface adsorption, thus representing a sig- 
nificant destabilization of the surface. Subsequent in- 
vestigations of this system in more detail indicated that 
structures involving on-surface and sub-surface oxygen 
functioned as nucleation sites and metastable precursor 
structures for bulk oxide formation and an atomic path- 
way for the transition to RUO2 was predicted.^" Similar 
studies for 0/Rh(lll)^^ have also recently been reported, 
as well as studies of the initial oxidation of alloy surfaces, 
i.e., NiAl(llO).^^ From a recent trend study of oxygen 
adsorption and incorporation at the basal planes of Ru, 
Rh, Pd, and Ag, it was found that the coverage at which 
the onset of oxygen occupation of sub-surface sites oc- 
curs, correlates closely with that predicted at which bulk 
oxide formation begins, and that this occurs at progres- 
sively lower coverages for the metals more to the right in 
the periodic table. ^'^ The energy cost for lattice distortion 
was found to play a key role, where it is greater for the 
the metals more to the left in the periodic table (i.e., Ru 
and Rh). With regard to pure free-electron- like metals, 
the 0/Al(lll)24 and O/Mg(0001)25 systems have been 
investigated as well by first-principles calculations. Ear- 
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lier studies using an embedded atom-type approach, fo- 
cussed on the energy barriers for O penetration through 
the surface layer of various metals and oxides, which is an 
important consideration particularly with regard to the 
kinetics.^® Our present investigation into the interaction 
of oxygen and Ag(lll), in addition to shading light on 
the understanding of how silver function as an efficient 
oxidation catalyst, also addresses the initial stages of ox- 
ide formation. Our general findings could have conse- 
quences for gold, also a noble metal oxidation catalyst^^ 
which, for the (111) surface, exhibits a restructuring and 
chcmisorption of oxygen atoms at elevated temperatures 
(500 - 800 K) and atmospheric pressure, ^^'^^ very simi- 
lar to Ag(lll). The findings of the present study could 
also have consequences for copper, which catalyses the 
oxidation of methanol to formaldehyde.^*^ 

Experimental evidence indicates that the (111) orien- 
tation is an important crystal face for real silver cata- 
lysts since at high temperatures, facets with this face re- 
sult -'^^^-'^^ due to it having the lowest surface energy. There 
are only two ordered phases of oxygen on Ag(lll) that 
have been reported; the (4 x Af^ and (\/3 x \/3)i?30° 
structures. '''^^ The latter actually exhibits a superstruc- 
ture given in matrix notation as (26 x 1;— 1 x 26). 
The (4 X 4) structure has recently been investigated by 
scanning tunneling microscopy (STM) where the atomic 
structure is proposed to involve a thin surface-oxide 
layer. ^^'^^ In this work (Ref. 33), it was furthermore pro- 
posed that at low coverage {9 = 0.05±0.03 ML), O atoms 
occupied sub-surface octahedral sites below the first Ag 
layer. As will bo shown in our present work, however, on- 
surface adsorption at low coverage is significantly more 
favorable. 

The (V3 X v^)i?30° phase has also been investigated 
by STM, as well as by reflection electron microscopy 
(REM) and reflection high-energy electron diffraction 
(RHEED). It forms only after exposure of silver (irre- 
spective of the crystal face) to O2 at atmospheric pressure 
and high temperatures (800-900 K).^'^^ Various atomic 
structures have been proposed for this phase; our recent 
STM simulation supported a surface-substitutional ge- 
ometry, however it was noted that this structure was only 
metastable?^ To date, the atomic geometry of neither of 
these ordered phases have been unambiguously confirmed 
which calls for further quantitative structural analyses. 

In our previous work,'^'* we studied both on-surface 
oxygen and oxygen adsorption in surface-substitutional 
sites on Ag(lll) for a wide range of coverages. The 
bondstrength of on-surface oxygen was found to be weak 
compared to the transition metals to the left of silver 
in the periodic table (e.g. Ru). With increasing cover- 
age, there is a strong repulsion between on-surface oxy- 
gen atoms, and on-surface oxygen adsorption becomes 
energetically unstable with respect to gas phase molec- 
ular oxygen for coverages greater than ~0.5 ML. In the 
present paper we investigate pure sub-surface oxygon and 
structures involving both oxygen adsorbed on the surface 



and in the sub-surface region, hereafter denoted as "on- 
surface+sub-surface" . We also investigate the proposed 
(4 X 4) surface-oxide phase, and oxide-like structures con- 
taining higher oxygen concentrations. These results, to- 
gether with our earlier study, provide a comprehensive 
understanding of the electronic and atomic structure, and 
stability of oxygen species at the Ag(lll) surface. 

The paper is organized as follows: The calculation 
method is briefly explained in Sec. 11, and in Sec. Ill, 
results are presented for purely sub-surface O structures 
(i.e. without on-surface oxygen) where the energetics, 
atomic geometry and electronic properties for different 
sites and for a wide range of coverages are described. 
Also energy barriers for oxygen penetration into the sub- 
surface region are reported. In Sec. IV structures are 
considered that involve both on-surface and sub-surface 
oxygen species, as well as the earlier proposed (4 x 4) 
geometry.^^ In Sec. V, geometries involving higher con- 
centrations of oxygen are considered, where the energetic 
preference of oxide-like structures is found which have a 
geometry similar to the (111) surface of Ag20. In this sec- 
tion, we also investigate an ozone-like species adsorbed 
at an Ag vacancy, as has been proposed to be the ac- 
tive species for ethylene epoxidation. The conclusions 
are given in Sec. VI. 

II. CALCULATION METHOD 

The density-functional theory (DFT) total energy 
calculations are performed using the pseudopotential 
plane-wave mothod'^'"^ within the generalized gradient ap- 
proximation (GGA) for the exchange-correlation func- 
tional. ^^'^^ The pseudopotentials are generated by the 
scheme of Troullier and Martins with the same func- 
tional. ^^'^^ The wave-functions are expanded in plane- 
waves with an energy cutoff of 50 Ry and the surface is 
modeled by a five layer slab separated by 15 A of vac- 
uum. Oxygen is placed on one side of the slab where the 
induced dipolc moment is taken into account by applying 
a dipolc correction.^" The positions of the oxygen atoms 
and the top three silver layers are relaxed until the forces 
on the atoms are less than 0.015 cV/A. In the (1 x 1) 
Ag(lll) surface unit cell, 21 special k-points are used 
in the surface irreducible Brillouin zone (IBZ) for the 
Brillouin-zone integration.'*^ These k-points are equiva- 
lent for all of the surface structures studied in order to 
maximize the accuracy when comparing the energies of 
different coverages as calculated in different supercells. 
We employ a Fermi function with a temperature broad- 
ening parameter of T'^' = 0.1 eV, and the total energy 
is extrapolated back to zero temperature. Results of de- 
tailed convergence tests can be found in Ref. 34. 

The average adsorption energy per oxygen atom, E^d, 
is calculated according to. 
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E,^ = _^[i;0/Ag(lll) _ (^Ag(lll) ^ ^o^free-O^] ^ 

No 

(1) 

where A^o is the total number of oxygen atoms of 
the adsorbatc-substrate system, and the total energy of 
the adsorbate-substrate system, the clean surface, and 
the free oxygen atom are represented by £;0/^s(m)j 
£:Ag(iii)^ and ^;frcc-o^ respectively. That is, i;0/Ag(iii) 
represents the 0/Ag system under investigation, which 
e.g. may involve purely sub-surface O, or on-surface and 
sub-surface species, or oxide-likc structures. 

To study the interaction between on-surface and sub- 
surface oxygen, we consider how the stability of on- 
surface O is affected by sub-surface O (and vice- versa). 
To do this, we define a so-called "removal energy" , which 
is the energy required to remove an on-surface O atom 
into the vacuum. It is given as, 

^removal _ _J^O/Ag(lll) _ ^^08ub/Ag(lll) ^ ^free-O^j 

(2) 

where ^03„b/Ag(iii) jg ^j^g ^^^3,1 energy of the "reference" 
system, i.e., that containing only (one or more) sub- 
surface O species, ijromovai ^^-^ .^^so bc thought of as the 
adsorption energy of an O atom onto the substrate which 
contains the sub-surface O atoms. An obvious analogous 
equation holds for the removal energy of a sub-surface O 
atom, 

_£;^emoval ^ _J^O/Ag(lll) _ ^^Oo„/Ag(lll) _|_ ^free-O-jj 

(3) 

where £;Oon/Ag(iii) jg ^j^g total energy of the "refer- 
ence" systcmi. i.e., that containing only the on-surface 
O species. We consider El^°^'^^ only for structures hav- 
ing one O atom on the surface and one below it. The 
adsorption and removal energies are defined such that a 
positive number indicates that the adsorption is exother- 
mic (stable) with respect to a free O atom, and a negative 
number, endothermic (imstable). 

The definitions of quantities that we use to analyze 
our results, e.g., total and projected density of states 
(PDOS), density difference distributions, and the work 
function, can be found in Refs. 42, 34 

III. PURE SUB-SURFACE OXYGEN 

We first study sub-surface oxygen without the presence 
of on-surface oxygen. For oxygen occupation in the sub- 
surface region there two possible sites: (i) the octahedral 
site, denoted hereafter as "octa" , which has six nearest- 
neighbour Ag atoms, three above and three below, and 
(ii) the tetrahedral site, which has four nearest-neighbour 
Ag atoms. Due to the presence of the surface, there are 
two types of tetrahedral sites; one is where there are three 
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FIG. 1. Calculated adsorption energy of oxygon (with re- 
spect to atomic oxygen) under the first Ag(lll) layer in octa 
(circles), tetra-I (diamonds), and tctra-II (triangles) sites, as 
a function of coverage. The horizontal full line represents half 
the experimental value of the binding energy of O2 (2.56 eV). 
The solid linos connecting the calculated adsorption energies 
are to guide the eye. 

silver atoms above it and one below, denoted as tetra- 

I, and the alternative one, tetra-II, is just the opposite 
with one surface Ag atom directly above, and three below 
it in the second Ag layer, as depicted in the insets of 
Fig. 1. We performed calculations for oxygen in these 
different sites for a wide range of coverages, i.e., 0.11 to 
1 ML. We focus on adsorption immediately below the 
first Ag layer as we find that oxygen adsorption deeper 
in the bulk is less favorable in every case (see Tab. I). In 
the following subsections, the energetics, the atomic and 
electronic structure, and energy barriers for O diffusion 
into the substrate are discussed, respectively. 

A. Energetics 

The calculated adsorption energies, i^adj for Ooctai 
Otetra-i, and Otetra-iii with respect to atomic oxygen, 
are plotted in Fig. 1 and are listed in Tab. I. For Oocta 
(circles), it can be seen that the adsorption energy de- 
creases when the coverage varies from 0.11 to 0.33 ML, in- 
dicating a repulsive interaction between the Oocta atoms; 
however, it increases from 0.33 to 0.50 ML, so that the 
effective interaction betwcien the Oocta atoms is attrac- 
tive. For coverages higher than 0.50 ML, the adsorption 
energy decreases again. Despite this dependence on cov- 
erage, we note that the overall variation in the magnitude 
of the adsorption energy is rather small, i.e. not greater 
than 0.3 eV - in sharp contrast to what we found for ad- 
sorption on the surface, where the corresponding value is 
1.93 eV.^'* Compared to Oocta, the tetrahedral sites are 
energetically less favorable for the whole coverage range. 
It can be noticed that the tetrahedral sites exhibit a sim- 
ilar dependence on coverage as Oocta- 
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TABLE I. Calculated energies (in eV) and first metal in- 
terlayer expansion for various coverages of sub-surface oxygen 

in the octa, tctra-I, and tctra-II sites. Sad is the adsorption 
energy with respect to atomic oxygen and i5b is the binding 
energy on the pre-distorted surface, and Ed is the energy cost 
of the distortion of the substrate as induced by sub-surface 
oxygen (without the presence of O, see text). Adi2 is the per- 
centage change in the first metal interlayer distance, where 
the center-of-mass of each layer is used. The calculated bulk 
interlayer distance is 2.43 A . Where available, values of the 
adsorption energy are given for O occupation under the sec- 
ond metal layer, i5ad(2nd). 
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From our calculations the relative stability of the var- 
ious systems are well defined, the absolute value of the 
adsorption energies is less so due to the error of stan- 
dard DFT calculations in describing the free O atom and 
molecule.^'* In particular, the experimental binding en- 
ergy of O2 per O atom is 2.56 eV while the calculated 
value is 3.16 eV/atom,^^ that is, the theoretical value is 
significantly overestimated. The deviation to the experi- 
mental value highlights the difficulty in obtaining quanti- 
tatively, the absolute adsorption energies. This is partic- 
ularly problematic for this system as can be appreciated 
from Fig. 1: On consideration of the theoretical value of 
3.16 eV, it would indicate that all of the sub-surface sites 
are unstable with respect to the oxygen molecule, while 
the experimental value indicates that the octa site is sta- 
ble with respect to O2 for all the coverage range. The 
adsorption energy of oxygen on the substrate is expected 
also to be overbound, although not nearly so severely as 
the O2 molecule. 

Compared to sub-surface adsorption, on-surface ad- 
sorption is notably more favorable for low coverages: 
compare 3.61 eV^* (fee) to 2.93 eV (octa. Tab. I) for 
coverage 0.11 ML and 3.52 eV^'' (fee) to 2.85 eV (octa, 
Tab. I) for coverage 0.25 ML. Due to the build-up of a 
strong repulsion between on-surface adsorbates,^^ for a 
coverage of 0.50 ML the energies of on-surface and sub- 
surface oxygen become practically degenerate (compare 
2.92 eV^^ (fee) and 2.94 eV (octa, Tab. I)). For higher 



coverages, on-surface adsorption becomes considerably 
less favorable than sub-surface oxygen 

The electron density at sub-surface sites is significantly 
greater than for on-surface sites which gives rise to a 
strong kinetic repulsion due to the orthogonalization of 
the wavefunctions of the substrate caused by the presence 
of oxygen. To reduce this repulsion, the lattice expands 
to an optimum value.^^ In order to aid our discussions 
concerning on-surface versus sub-surface adsorption, it is 
constructive to regard the sub-surface adsorption energy, 
i^ad) as consisting of two contributions: (i) the energy 
cost of distorting the (clean) substrate lattice to the ge- 
ometry induced by sub-surface oxygen, E^, and (ii) the 
binding energy of oxygen on the pre-distorted substrate, 
E\), i.e., = E\) — Ed- We may then regard the bind- 
ing energy as being due only to electronic effects, as we 
have removed the energy cost of distorting the lattice; 
that is, as purely reflecting the bond strength of the O- 
Ag interaction. We have calculated these contributions 
for the octa and tetra-I sites which are given in Tab. I 
and are also plotted in Fig. 2, along with the correspond- 
ing on-surface adsorption energies^^ for comparison. In 
addition, the percentage change in the first metal-metal 
interlayer distance is shown. 

Comparing the binding energies of O on the pre- 
distorted surface with on-surface adsorption (see Fig. 2a), 
it can be seen that for coverages 0.25 and 0.33 ML, on- 
surface adsorption is still (slightly) energetically more fa- 
vorable than sub-surface adsorption in either the octa or 
tetra-I sites, thus the preference for O to be on the sur- 
face is not due to the energy cost of distortion of the 
substrate for the case of sub-surface adsorption. It can 
be seen, however (from e.g. Fig. 2b) that the energy 
cost of distortion is significant and does reduce the ad- 
sorption energy of sub-surface oxygen substantially. For 
example, at 0.50 ML the binding energy of sub-surface 
oxygen in the octa site on the pre-distorted substrate 
is 0.48 eV greater than adsorption on the surface, but 
the adsorption energy (i.e. when paying for the distor- 
tion cost) is practically the same. Furthermore, it can be 
seen that like the adsorption energy, only to a greater de- 
gree, the binding energy of sub-surface oxygen decreases 
in the coverage range 0.50 - 1.0 ML. This repulsive inter- 
action is attributed to unfavorable 0-0 interactions due 
to the close proximity of the partially negatively charged 
O atoms to one another and the larger expansion of the 
metal-metal spacing for these coverages, which less effec- 
tively screens the O atoms. 

The energy cost of distorting the substrate (Fig. 2b) 
per O atom generally decreases with increasing O cov- 
erage due to an effective sharing of the energy cost by 
neighboring O atoms. The exception, from coverage 
0.33 to 0.50 ML, where the distortion energy slightly 
increases, is because there is a significant increase of 
the metal interlayer spacing for the latter coverage (see 
Fig. 2c). The distortion energy is always larger for the 
tetra-I site, which is understandable due to the larger 
oxygen-induced expansion of the first metal interlayer 
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FIG. 2. (a) Binding energies of O on the pre-distorted sub- 
strate, Eb, in the octa (filled circles) and tetra-I sites (filled 
squares) and the corresponding adsorption energies, E^d, of O 
in the octa (open circles) and tctra-I (open squares) sites. The 
adsorption energies of on-surface fcc-O are also shown (open 
diamonds), (b) Energy cost, E^, of distorting the substrate 
lattice to the geometry induced by sub-surface oxygen (but 
without the presence of oxygen) for the octa (open circles) 
and tetra-I (open squares) sites, (c) Percentage change of the 
first metal interlayer distance compared to the bulk value for 
the octa (circles) and tetra-I (squares) sites. 




FIG. 3. The difference electron density, n(r)^, (see the def- 
inition in Ref. 34) for oxygen in the sub-surface tetra-I site 
for coverage 0.25 ML. The contour plane is in the [211] direc- 
tion and is perpendicular to the Ag(lll) surface. The unit 
is 10^^ Bohr^^. The positions of the O and Ag atoms are 
indicated by the arrows. 



due to a less efficient sereening by the metal atoms of the 
partially negatively charged sub-surface O atoms, due to 
the larger Ag-Ag interlayer spacing of the tetra-I site as 
shown clearly in Fig. 2c. 

The above findings highlight the interplay between the 
energy gain due to bond formation and the energy cost 
due to structural deformation. These effects are mate- 
rials dependent: To the left of Ag in the periodic table, 
the transition metals are "harder" (i.e., have higher bulk 
moduli), and the energy cost of deformation is higher;^^ 
on the other hand, the energy gain due to bond forma- 
tion between oxygen and the transition metal is stronger 
as well. 



spacing (see Tab. I and Fig. 2c) caused by the smaller 
space in which to accommodate the sub-surface O atom. 

As seen clearly from Fig. 1, the adsorption energy for 
the octa site is greater than the tctrahcdral sites for all 
the coverages considered. If we remove, however, the 
distortion energy, we see from Tab. I (and Fig. 2a) that 
for coverage 0.25 ML, the tetra-I site has in fact slightly 
greater binding energy (compare 3.49 to 3.40 eV). In this 
case it is the larger cost of distortion for the tetra-I site 
(compare 0.88 to 0.54 eV) that results in the octa site 
having the larger adsorption energy. For all the higher 
coverages, however, the octa site has the stronger 0-Ag 
binding and even though the energy cost of distortion 
for the tetra-I site is larger than the octa site in all cases, 
this does not determine the adsorption site preference 
and make the tetra-I site less favorable than the octa site. 
The reason that the binding energy of the tetra-I site is 
less favorable than the octa site for coverages greater than 
0.25 ML, is therefore an electronic effect. This could be 



B. Atomic structure and electronic properties 

As seen above, sub-surface oxygen occupation induces 
a significant structural distortion of the Ag(lll) sub- 
strate as shown by the variation of the first interlayer 
spacing in Tab. I and in Fig. 2c: Namely, the expansion 
varies from 0.16 to 0.81 A for Oocta when 9 increases from 
0.11 to 1.00 ML, and from 0.38 to 1.14 A for Otetra-i, 
and from 0.43 to 1.20 A for Otetra-ii when 6 increases 
from 0.25 to 1.00 ML. This is in sharp contrast to what 
we found in our study of on-surface oxygen adsorption 
where the contraction of the first interlayer spacing in- 
duced by on-surface oxygen was less than 0.02 A for the 
whole coverage range investigated (0.11 to 1.00 ML).^'* 
The large expansion is, however, similar to what has been 
found for sub-surface O at other metal surfaces. 

The 0-Ag bondlengths are generally longer for the 
case of the octa site compared to the tetra sites, which 
is in-line with the understanding that adsorbates in 
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FIG. 4. Total projected density of states (PDOS) for oxy- 
gen in the sub-surface tetra-I site for coverage 0.25 ML. The 
Ag PDOS of the first and second metal layer atoms that are 
bonded to oxygon arc indicated by dashed and dot-dashed 
lines, respectively, and the O PDOS is indicated by a solid 
line. The vertical dashed line marks the Fermi energy. 



higher coordinated sites have longer bondlengths com- 
pared to lower coordination sites. "'^'"'^ Most of the 0-Ag 
bondlengths involving sub-surface oxygen are longer than 
the typical bondlength between on-surface oxygen and 
Ag (2.17 A), i.e., roughly 8.8, 2.3, 2.8% longer for the 
octa, tetra-I, and tetra-II sites, respectively. 

The change in the work function induced by sub- 
surface oxygen is modest, i.e., less than 0.56 eV for all 
of the adsorption sites, throughout the range of cover- 
ages considered, which is in strong contrast to the case 
of on-surface oxygen adsorption, where it is as large as 
4 cV at 1.00 ML.34 This is due to the location of O un- 
der the surface Ag layer, which screens well the partially 
negatively charged oxygen atoms, and in addition there 
is, to a large degree, a cancellation of dipole moments 
involving oxygen and the top two Ag layers as indicated 
by the difference electron density in Fig. 3. 

The electron density of the exposed (surface) silver 
atoms and the silver atoms bonded to O in the second 
layer, is depleted, while there is a large accumulation on 
the O atom. The result is very similar for the octa site. 

Similarly to on-surface oxygen adsorption,^^ the inter- 
action between sub-surface oxygen and silver is mainly 
via hybridization of the 0-2p and Ag-Ad-bsp states, where 
antibonding states are largely occupied, indicating an 
ionic nature and explaining the relatively weak binding 
energy compared to the metals to the left of Ag in the 
periodic table. This can be seen from the total projected 
density of states of the O atom and of the Ag atoms 
in the first and second layers that are bonded to O, as 
shown in Fig. 4 for oxygen in the tetra-I site at coverage 
0.25 ML. The band width of the Ag-M states are no- 
tably narrower for the surface Ag atoms due to the lower 
coordination. As for on-surface oxygen adsorption, no 
other peaks occur between the 0-2s state and the lower 



edge of the Ag-4d band; we point this out because such 
a feature has been found in some experimental studies, 
the origin of which is under debate (see e.g., Ref. 2). 
Clearly, purely sub-surface oxygen is not responsible for 
these features, nor is purely on- surface oxygen or surface- 
substitutional oxygen.^** This will be discussed in more 
detail in Sec. VB. Compared to on-surface oxygen, the 
0-2s PDOS for sub-surface oxygen is at a lower energy, 
which, from the initial state theory of core-level shifts, in- 
dicates that sub-surface oxygen is less negatively charged. 



C. Energy barriers for O penetration into the 
substrate 

To investigate sub-surface oxygen in more detail with 
regard to the energetics of its formation, we study the 
penetration of O through the surface Ag layer from the 
on-surface (fee) site to the sub-surface octahedral site at 
low coverage (0.11 ML). The transition state corresponds 
to the oxygen atom in the plane of the surface Ag(lll) 
layer. We fully relax the first three Ag layers, but to stay 
at the transition state, we restrain O to be in the plane of 
the first Ag(lll) layer (if we did not do this, the O atom 
would move to one of the two minima either side). The 
vertical distance of this top layer, however, is relaxed. 
We also considered the diffusion barriers for oxygen un- 
der the first Ag layer to move deeper into the bulk, i.e. 
to diffuse through the second Ag layer by analogous cal- 
culations. In this case the pathway is from the octa site 
under the first Ag layer to the tetra-I site under the sec- 
ond Ag layer. Oxygen will then move to the energetically 
more favorable octa site under the second Ag layer. We 
note that due to the use of a (3 x 3) surface cell, which 
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FIG. 5. Adsorption energy of oxygen at a coverage of 
0.11 ML for pcuotratiou into the sub-surface octa site from the 
on-surface fcc-hoUow site through the first Ag layer, and from 
the octa site to the "bulk" tetrahedral site, through the sec- 
ond Ag layer (from right to left). From the tetra-I site O will 
move to the more favorable "bulk" octa site under the second 
Ag layer. The short lines indicate the adsorption energy with 
respect to half of the binding energy of the oxygen molecule 
(experimental value). "TS" represents the transition states, 
and the numbers on the arrows indicate the energy barriers. 
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will restrict complete lateral relaxations at the transition 
state, the barriers could be slightly lower. The present 
values may be taken as an upper limit. Prom Fig. 5 it can 
be seen that the energy barrier for fcc-oxygen to enter the 
sub-surface region is 0.86 eV. For the reverse process, i.e., 
to diffuse to the surface from the octa site, the barrier 
is only 0.18 eV. Assuming that the prefactor for the two 
diffusion directions is the same, and using a simple Ar- 
rhcnius equation, the ratio of the diffusion coefficients for 
movement from the surface to the sub-surface, compared 
to from the sub-surface to the surface is: 5.1x10"^ at 
470 K and LGxlO"" at 900 K. These kinetic considera- 
tions indicate that chemisorbed oxygen will stay on the 
surface at low coverage and at temperatures of around 
470 K. Even at the higher temperature of 900 K, the 
difference is still significant. 

The obtained energy barrier to move deeper into the 
bulk from under the first Ag layer to under the second Ag 
layer is 0.92 eV. This value is slightly larger compared to 
the case of oxygen diffusing through the first metal layer; 
this is due to the presence of the surface in that it costs 
less energy for the surface Ag atoms to relax compared 
to the Ag atoms in the second layer, which are more 
constricted due to other Ag atoms both above and below. 
Although configurational entropy for O in the bulk could 
drive oxygen to diffuse deeper into bulk region, the above 
results clearly show that there is large potential gradient 
indicating that the oxygen, in the lower coverage regime, 
will be concentrated near surface. 

As noted in the introduction, a chemisorbed oxygen 
species at coverage 0.05±0.03 ML was observed, by scan- 
ning tunneling microscopy (STM).^'^ The associated fea- 
ture exhibited a depression with a diameter of 9±lA 
which, according to the STM stimulation, suggested that 
the oxygen atom is located in the octa site under the first 
silver layer. On-surface adsorption was excluded since 
the diameter of the depression in the stimulations was 
significantly less, i.e. just 4±lA. The tetra-I site was ex- 
cluded on the basis that it offered less space and was ex- 
pected to be energetically unfavorable. The alternative 
tetra-II site was not mentioned. No calculations, how- 
ever, of energies were performed in this study. As our 
results show, at coverage 0.11 ML, the binding energy 
of O on the surface is energetically significantly more fa- 
vorable than sub-surface oxygen which docs not support 
this assignment. The effective 0-0 repulsion is not very 
strong in the coverage range 0.11 to 0.25, i.e. the adsorp- 
tion energy is only 0.09 eV smaller for the latter coverage. 
Thus we believe that for the lower coverage of 0.05 ML, 
the energetic preference for on-surface adsorption will not 
change. 

IV. EFFECT OF SUB-SURFACE OXYGEN 

In the previous section it was shown how on-surface 
oxygen was energetically favorable for coverages less than 



w0.50 ML compared to pure sub-surface adsorption, but 
for higher coverages, sub-surface adsorption was pre- 
ferred. It was also seen how the presence of sub-surface 
oxygen changes the nature of the surface Ag atoms, 
namely, by causing a depletion of the electron density. 
In order to investigate the effect that sub-surface oxygen 
has on the reactivity of the surface, as a first step we 
study the adsorption of additional oxygen on the surface 
with sub-surface O present. This study is interesting in 
the sense that if on-surface oxygen is taken to represent 
an arbitrary electronegative adparticle, its behavior may 
reflect a general tendency. Investigating the interaction 
between on- and sub-surface oxygen is also essential in 
relation to understanding the initial formation of surface 
oxides. We first consider many different atomic arrange- 
ments involving structures with an O atom on the surface 
and one under the surface Ag layer in order to identify 
energetically favorable geometries. 

A. Energetics 

We consider the case of one on-surface oxygen atom 
in the fee site, Ofcc, of the (2 x 2) cell plus one sub- 
surface oxygen atom, i.e., a total coverage of 0.50 ML. 
With this arrangement, there are three high- symmetry 
sub-surface sites, namely, octa, tetra-I, and tetra-II, as 
described above, but due to the presence of Ofcc, there are 
now two different kinds of each of these sites, which can 
be differentiated by their relative distance to Ofcc. We 
therefore assign the label "Inn" (1st nearest-neighbor) 
to indicate the three sub-surface sites closest to Ofcc and 
"2nn" to indicate those that are the second nearest neigh- 
bor sub-surface sites. These six possible sub-surface sites 
are shown in Fig. 6. We tested all of these possibilities, 
where the positions of the oxygen atoms and top three 
silver metal layers were fully relaxed. Similar geometries 
exist when on-surface oxygen is located in the hep-hollow 
site; in the light of what is learnt from the results for O 
in the fee site however, we viewed it necessary to only 
consider two structures for this case. 

We first discuss results for sub-surface adsorption with 
on-surface O in fee sites. The calculated (average) ad- 
sorption energies are given in Tab. II, along with the 
removal energies (cf. Eqs. 2 and 3), and the 0-Ag 
bondlengths and metal interlayer expansions, as well as 
the work function change. 

It can be seen that the values depend strongly on the 
location of sub-surface oxygen: The average adsorption 
energies vary from 2.77 to 3.27 eV. Furthermore, the Inn 
set of structures are energetically unfavorable compared 
to the 2nn set. This is due to the close proximity of the 
former to Ofcc which causes a considerable 0-0 repulsion. 
For example, the adsorption energy of the most favorable 
structure of the Inn set is 2.92 eV, while the least favor- 
able structure of the 2nn set is 2.99 eV. The energetically 
most favorable structure of all is that where sub-surface 
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fcc/octa (Inn) fcc/tetra-I (Inn) fcc/tetra-II (Inn) 




fcc/octa (2nn) fcc/tetra-I (2nn) fcc/tetra-II (2nn) 
FIG. 6. Schematic geometry of the (2 x 2) surface unit 
cells with on-surface oxygen in the fee site, Ofcc , (larger dark 
circles) for a coverage of 0.25 ML. The various sub-surface 
sites are depicted as the small dark circles and the resulting 
on-surface+sub-surface geometries defined by the label (e.g. 
fcc/octa (Inn)). There are two different types of sub-surface 
sites; ones that are closest to Ofcc, denoted 'Inn' and ones that 
are further away, denoted '2nn'. Note, as indicated by the 'x' 
on the atoms, the Inn-Oocta site is located directly below 
Ofcc and the tetra/-II site is located directly below surface Ag 
atoms. 

oxygen is located in the tetra-I site of the 2nn set. This 
structure is denoted as (Ofcc/Otetra-i)0=o.5O and can be 
seen more clearly in Fig. 10b. 

To investigate the affect of sub-surface oxygen on the 
bond strength of Ofcc, wc consider the removal energy 
cf. Eq. 2. We calculated this quantity for O in all of the 
various sub-surface sites. For sub-surface O in the octa, 
tctra-I, and tetra-II sites of the Inn set, the respective re- 
moval energies are: 2.98, 2.94 and 3.36 eV. For the same 
sites but in the 2nn set, the values are 3.10, 3.93 and 
3.64 eV, respectively. These vahies can be compared to 
the adsorption energy of Ofcc (at 0.25 ML) without the 
presence of sub-surface oxygen which is 3.52 eV.^^ Thus, 
depending on the location of sub-surface oxygen, it can 
either destabilize or stabilize significantly the interaction 
between on-surface oxygen and the metal substrate. The 
energetically most stable structure, (Ofcc/Otetra-i)6i=o.50) 
stabilizes on-surface O by 0.41 eV (compare 3.93 to 
3.52 cV). The affect of the presence of on-surface O on the 
binding of sub-surface O (cf. Eq. 3) is like that found for 
the affect of sub-surface O on the binding of on-surface O. 
This can be seen from inspection of the values of £"3""°™' 
in Tab. II, together with the values of i?ad in Tab. I for 
structures involving purely sub-surface oxygen. 

For on-surface oxygen in the hcp-hollow site, we 
considered the octa and the tetra-I sites of the 
2nn set for sub-surface O. The former, denoted as 
(Ohcp/Oocta)6i=o.50, IS geometrically very similar to the 
energetically favorable structure, (Ofcc/Otetra-i)0=o.5O, 
namely, both involve a linear 0-Ag-O bonding arrange- 
ment. And as can be seen from the results presented 
in Tab. II, it is also a low energy structure; in fact, it 
is the lowest of all. That it is slightly more favorable 



TABLE II. Calculated atomic geometries, energetics and 
work function change for oxygen adsorption in structures in- 
volving both on-surface fcc-hoUow sites and sub-surface sites 
at a total coverage of 0.50 ML. is the bondlength be- 
tween on-surface oxygen and silver atoms in the first metal 
layer, while Rf^^ is the bondlengths between sub-surface oxy- 
gen and silver atoms in the first metal layers. di2 is the first 
metal interlayer distance given with respect to the center of 
gravity of each layer. The bulk interlayer distance is 2.43 A . 
£'ad is the average adsorption energy with respect to atomic 
oxygen, while and EIX°""^ are the removal energies 

for on-surface and sub-surface oxygen atoms, as defined in 
Eqs. 2 and 3. A<I> is the change in work function with respect 
to the clc^an surface (cak-ulatcd to bo 4.45 oV). The unit of 
length is A and energy is e\'. 
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1.02 
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2.16 


2.18 


2.86 


2.77 


2.94 


2.14 


1.59 



(by 0.06 eV) may be understood due to the fact that the 
pure sub-surface Oocta structure at 0.25 ML is energeti- 
cally more favorable than the Otctra-i site (by 0.24 cV) 
as shown in Tab. I, while the difference between pure on- 
surface Ofcc and Ohcp is less (0.11 eV). 

Comparing the average adsorption energies of the 
structures described above for coverage 0.50 ML to the 
structures at same coverage, but where oxygen occupies 
exclusively either on-surface (£ad = 2.92 eV for fee oxy- 
gen) or sub-surface sites (E'ad = 2.94 eV for Oocta); it 
is found that the mixed structures, (Ofcc/Otctra-i)6i=o.50 
(3.27 eV) and (Ohcp/Oocta)e=o.50 (3.33 eV), are signifi- 
cantly energetically more favorable. The average adsorp- 
tion energies of these mixed structures are, however, not 
greater than that of on-surface adsorption at the lower 
coverage of 0.25 ML. Other geometries are also possible 
for coverage 0.50 ML, for example, pure on-surface O in- 
volving hep and fee sites, (Ofcc + Ohcp)e=o.50, and pure 
sub-surface O in different sites, (Oocta + Otetra-ii)e=o.50 
and (Otetra-i + Otetra-ii)e=o.50, where the two sub- 
surface oxygen atoms are separated as far as possible 
from each other. All these geometries were tested and 
found to be energetically unfavorable, as shown by their 
respective average adsorption energies, namely, 2.41, 2.94 
and 2.95 eV (cf. Tab. III). Altogether, our calculations 
indicate that oxygen will enter the sub-surface region for 
coverages greater than 0.25 ML. We note that this cover- 
age could be less, but to determine this would require the 
use of larger supercells; the value of ^ = 0.25 ML may 
therefore be regarded as an upper limit. Clearly this is 
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in contrast to O adsorption on Ru(OOOl) where O ad- 
sorption in sub-surface sites only becomes energetically 
preferred after completion of a full monolayer. -"^^'-^^ 

Viewing oxygen as a typical electronegative atom, its 
function at the sub-surface site may be replaced by other 
electronegative species, like chlorine, which is actually 
used as a promoter in some reactions that silver catal- 
yses.^ It may be expected then that chlorine will also 
stabilize on-surface electronegative species. Also, the 
found enhanced stability of on-surface oxygen due to sub- 
surface O may also reflect the behavior of other similar 
species at the surface. This could explain the reported 
stabilization of CO on Ag(lOO) which was attributed to 
sub-surface oxygen.^^ 

B. Atomic structure and electronic properties 

In the following we describe the main structural and 
electronic properties of the energetically most favorable 
on-surface-|-sub-surface structures. Full details concern- 
ing the 0-Ag bondlengths and metal interlayer spacing 
arc; given in Tab. II. The most notable feature of the low 
energy structures is a characteristic linear Oon-Ag-Ogub 
bonding arrangement of the uppermost atoms. In-line 
with their energetic preference, they also have shorter 
0-Ag bondlengths; for example, for (Ofcc/Otetra-i)e=o.50 
the Ofcc-Ag bondlength, R°'^, and the bondlength of sub- 
surface O to the first Ag layer, Rf"^, are 2.10 and 2.09 A , 
respectively, and for (Ohcp/Oocta)e=o.50j the analogous 
values are 2.15 and 2.11 A, respectively. These can be 
seen from inspection of Tab. II to be shorter than the less 
favorable on-surface-fsub-surface geometries. They are 
also shorter compared to the corresponding bondlengths 
of the pure on-surface and pure sub-surface structures, 
e.g., 2.17 A for on-surface fee oxygen, 2.12 for Otetra-i, 
and 2.24 A for Oocta at coverage 0.25 ML. Compared 
to the values of Rf^^ (bondlengths of sub-surface O to 
the top Ag layer), for the two most favorable structures, 
the bondlengths of sub-surface oxygen to the second Ag 
layer (not shown here) are pronouncedly longer, which 
indicates a stronger bonding of sub-surface O to the sur- 
face Ag layer, which is also bonded to a surface O atom. 

In Tab. II the work function changes with respect to 
the clean Ag(lll) surface due to O adsorption for the 
on-surface-|-sub-surface structures described above are 
listed. The values vary from 0.84 to 1.81 eV, and for 
the energetically most favorable structures the values are, 

(Ofcc/Otetra-l)e=0.50 1-15 eV and (Ohcp/Oocta)e=0.50, 

1.02 eV, respectively. For comparison, the change in work 

function for on-surface oxygen at coverage 0.25 ML are 
1.23 eV (fee) and 1.31 eV (hcp),*^^ and the values of the 
work function for the on-surface fee and hep sites at the 
same coverage of 0.50 ML are more than twice as large 
as the structures involving both on- and sub-surface oxy- 
gen. Sub-surface O competes with on-surface oxygen 
for the bonding charge of the surface Ag atoms, as a con- 



sequence, the on-surface O species are slightly less nega- 
tively charged compared to when there is no sub-surface 
O. This can seen by comparing Fig. 8 (below) with that 
of Fig. 4 in Ref. 34, and is also indicated by the lower 
energy position of the 0-2s levels for the favorable on- 
surface-|-sub-surface structure as discussed below. 

In Fig. 7 we show the projected density of states for 

Ofcc and Otetra-I 

atoms of the energetically favorable (Ofcc/Ototra-i)e=o.50 
structure. They are labeled "Ofcc(Ofcc/Otetra-i)0=o.5o" 
and "Otetra-i(Ofcc/Otetra-i)e=o.50 respectively, and are 
denoted by thick continuous lines in the upper and lower 
parts of the figure. We also show the PDOS for the 
first layer silver atoms which are bonded to these two 
O atoms (dotted lines in Fig. 7). For comparison we 
show the PDOS for the Ofcc and Oocta atoms of the ener- 
getically unfavorable (Ofcc/Oocta)e=o.50 structure (thick 
dashed curves in the upper and lower parts of the fig- 
ure, respectively), as well as the PDOS for the structures 
containing purely Ofcc (i-c, without Otetra-i) and purely 
Otetra-i (i-e., without Ofcc) atoms. 

Considering first the results for (Ofcc/Otetra-i)e=o.50, 
it can be seen that the PDOS for Ofcc of this structure is 
shifted to a lower energy compared to the pure on-surface 
Ofcc species (compare the thick and thin continuous lines 
respectively, in the upper plots). The antibonding states 
of Ofcc of the (Ofcc/0totra-i)e=o.5o structure are around 

0. 74 eV lower compared to the case of pure Ofcc- A sim- 
ilar shift occurs for the 0-2s state. For Otetra-i of the 
(Ofcc/Otctra-i)e=o.50 structure, as compared to Otetra-i 
of the pure sub-surface Otctra-i at 0.25 ML, it exhibits 
a similar behavior but the magnitude of the shift is less, 

1. e., around 0.56 eV for the valence PDOS and 0.45 eV 
for the 2s-state. Due to the mentioned shifts to lower 
energies, the DOS at the Fermi energy decrease for the 
(Ofcc/Otetra-i)e=o.50 systcm compared to that not con- 
taining sub-surface oxygen. This is also the case for 
the total DOS. The density of states for the other en- 
ergetically favorable structure, (Ohcp/Oocta)e=o.50) (not 
shown) are very similar to that of (Ofcc/Otetra-i)6i=o.50 
which is in-line with their similar energetics. 

We now consider an energetically unfavorable con- 
figuration for comparison, namely, (Ofcc/Oocta)0=o.5O, 
(where Oocta is in the 2nn set of Ofcc, see Fig. 6 and 
Tab. II). A different behavior to that described above is 
found: Here the PDOS for on-surface oxygen (Ofcc) shifts 
towards higher energies with the presence of sub-surface 
O, as does the 0-2s state (compare the thick dashed and 
thin continuous lines in the upper plot of Fig. 7). Fur- 
thermore, we find that all the energetically unfavorable 
structures exhibit a higher total DOS at the Fermi level 
compared to the energetically favorable ones. 

To gain further insight into the bonding mechanism, 
we calculate the difference electron density distribution 
for (Ofcc/Otetra-i)e=o.50, as shown in Fig. 8. A depletion 
of Ag-4(ixz.yz states of the surface silver atoms bonded to 
the Ofcc and Otetra-i oxygen atoms occurs, while there 
is an increase in electron density on the O atoms. The 



9 



> 



o 



a 

o 

<u 



5 

4 

3 

2 

1 


5 



3 - 











— o,„(o,^yo,„ 


ra-I^e=ll.5l) 


- 




^Icc, 6=0.25 








- ■ 0,„(0 J0„ 


ta)e-0.50 






Ag(OJO„ 


tra-l)e=0.50 


- 






Am I 




- 


^ Aa 




T; 




- 






/\A^ 1 




- 










- 




— o, 


„-.(o,„/o 


etra-l)e=0.50 






— 0, 


Ira-I, e=0.25 








- - o_ 






_ 




Ag(0,„/0„„ 




- 


1 '1 \' 










(1 1 \ 










li 












/ 'H 










' /''/■• 


W A\ 1 












"s. _j --- 



-8 



-6 



-4 



-2 







E-Ep (eV) 

FIG. 7. Total projected density of states (PDOS) 
for the favorable (0fcc/0tctra-i)e=().5() and unfavorable 
(Ofcc/Oocta)e=o.50 structures. PDOS for Ofcc (upper figure) 
and Otetra-i (lower figure) are shown by thick solid lines, and 
the silver atoms shared by the two oxygen atoms are shown 
as dotted lines. For comparison, Otcc without the presence of 
sub-surface oxygen (thin solid line) and with sub-surface oxy- 
gen in the octa site (thick dashed line) are shown in the upper 
figure. Also for comparison (lower figure) are the PDOS for 
Otetra-i without on-surfacc oxygen (shown as thin solid line) 
as well as sub-surface O in the octa site of the unfavorable 
structure (thick dashed line). The energy zero is the Fermi 
energy. 



coupling between Ototra-i and the silver atom directly 
below in the second metal layer is comparably weak, as 
indicated by the modest perturbation of electron density, 
and as also suggested by the longer bondlcngth described 
above. It is the Ofcc-Ag-Oictra-i linear coordination that 
allows a strong hybridization of the Ag-4(ixz,yz orbital 
with the Ofcc and Otetra-i species, and also screens well 
the electrostatic interaction between them. This is also 
the case for the (Ohcp/Oocta)6i=o.50 system (not shown) 
which contains the same linear coordination (i.e., Ohcp- 
Ag-Oocta) and weaker coupling to the substrate under- 
neath. 

The energetic preference for the fcc/tetra-I and 
hcp/octa (i.e. a linear 0-metal-O) arrangement is found 
also at higher coverage. In particular, we also considered 
oxygen coverages of 2.00 ML. These striictures involve a 
full monolayer on the surface and full monolayer under 
the first Ag layer, where there are a total of six possible 
configurations. We tested all of these and as mentioned 




FIG. 8. The difference electron density, n{r)^ (see the defi- 
nition in Ref. 34) for (Ofcc/Otetra-i)e=o.50. The contour plane 
is in the [211] direction and is perpendicular to the Ag(lll) 
surface. The unit is 10~^ Bohr~^. The positions of the O and 
Ag atoms are indicated by the arrows. 



above, the energetically favorable structure is again that 
involving on-surface O in the hep site and sub-surface O 
in the octa site,(0hcp/0octa)e=2.oo (by 0.07 eV), and the 
next most favorable structure is that for on-surface O in 
the fee sites and sub-surface O atoms in the tetra-I sites, 
(Ofcc/0iotra-i)e=2.oo- The average adsorption energies 
of these structures, are however, unstable with respect 
to gas phase O2. Such linear 0-metal-O geometries have 
also been found to be the most energetically favorable of 
those considered for other transition metals, e.g. Ru^° 
and Rh.^^ As will be seen below, this is actually the pri- 
mary building block involved for the energetically favor- 
able structures involving higher oxygen concentrations as 
well as in silver(I) oxide, Ag20. In addition, oxygen ad- 
sorption on Ag(llO) with (nx 1) reconstructions'^® involve 
a similar 0-Ag-O configuration. 




FIG. 9. Sketch of the atomic structure proposed in the 
STM study of Ref. 32 for the (4 x 4) phase. The surface unit 
cell is indicated. The oxygen atoms are represented by small 
dark circles, the uppermost Ag atoms by grey circles, and 
the intact plane of Ag(lll) atoms lying below the 'O-Ag-O' 
trilayer are represented as the open circles. 
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C. The (4 X 4)-0/Ag(lll) phase 

The proposed atomic geometry for the (4 x 4)- 

0/Ag(lll) phasc"^^''^'' bears a very close resemblance to 
the identified low energy (Ofcc/0tetra-i)e=o.50 structure. 
The difference being that the surface Ag atoms of the for- 
mer structure are not commensurate with the underlying 
(1 X 1) surface unit cell, i.e. the upper 0-Ag-O 'tri-layer' 
is laterally expanded to equal that of the (111) surface 
of bulk Ag20. In addition, to avoid Ag atoms sitting 
in top sites directly on top of the underlying Ag atoms, 
these Ag atoms are missing, i.e. one per (4 x 4) cell. 
This results in a stoichiometry of the tri-layer of Agi.gsO 
with a corresponding oxygen coverage of 0.375 ML (see 
Fig. 9). We calculated the average O adsorption energy 
of this structure (including full atomic optimization) and 
find that it has a low energy (see Fig. 14). In obtaining 
this value, we have used the chemical potential of Ag in 
bulk to take into account the energy of the Ag atoms in 
the mixed tri-layer. We also considered a 'stoichiometric 
structure' i.e., corresponding to Ag20 within the tri-layer 
in which the Ag atoms sitting on top of underlying sub- 
strate atoms are present and found, in agreement with 
other recent calculations of these two geometries,^^ that 
it is less favorable by 0.7 eV per cell. 



V. FORMATION OF A SURFACE OXIDE 
A. Accumulation of oxygen in the sub-surface region 

In order to gain insight into the stability of structures 
with higher concentrations of oxygen, we performed cal- 
culations for many different geometries. Since there are 
no experimental results to provide guidance, and no addi- 
tional ordered phases observed, we use (2x2) surface cells 
and systematically investigate the many possible atomic 
configurations at the given coverages. Such a study would 
be unfeasible with the larger (4 x 4) cells, and given the 
very similar local atomic geometry mentioned above of 
the proposed (4 x 4) structure and our identified low en- 
ergy structure at 0.5 ML (and similar electronic struc- 
ture), we view this as an appropriate strategy. 

We first investigate various geometries involving 
0.75 ML. Considering the addition of another O atom 
to the energetically favorable systems with total O cov- 
erage 0.50 ML, our results show that oxygen will oc- 
cupy sub-surface sites rather than on-surface sites to 
avoid the strong dipole-dipole repulsion. As seen in 
Tab. Ill, for the energetically favorable configuration de- 
scribed above for 0.50 ML of oxygen, (Ohcp/Oocta)e=o.50, 
the addition of an O atom to the on-surface hep site 
yields an average binding energy of 2.97 eV, but it is 
3.25 eV if the additional O atom occupies the sub-surface 
octa site below the first Ag layer, (Ohcp/Oocta)e=o.75- 
A similar result is found for the other energetically fa- 
vorable structure at 0.50 ML, (Ofcc/Otetra-i)e=o.50) in 



that it is energetically favorable for the additional oxy- 
gen atom to occupy a sub-surface site, namely the tetra- 
I site, (Ofcc/Otctra-i)e=o.75, (3.06 eV), instead of an 
on-surface fcc-hollow site (2.90 eV). The reason that 
the average binding energy of (Ohcp/Oocta)e=o.75 with 
0.50 ML octa oxygen is 0.19 eV more favorable than 
(Ofcc/Otetra-i)e=o.75 with 0.50 ML tctra-I oxygen, is 
mainly due to the difference in energy of the pure sub- 
surface Oocta and the pure sub-surface Otetra-i structures 



% ^ % ^ ^ 
(a) 3.52 eV/ 0.25 ML 



- 0.8 % 



0.0% 




19.0 % 



1.7 % 



(b) 3.27 eV/ 0.50 ML 




36.9' 



■ 3.4 % 



(c) 3.29 eV / 0.75 ML 




40.9 % 



37.6' 



(d)3.22 eV/ 1.25 ML 

FIG. 10. Atomic geometry of energetically favorable struc- 
tures for increasing oxygen concentrations: (a) on-surface fee 
oxygen, (b) fee oxygen plus sub-surface oxygen in the tetra-I 
site, (c) as for (b) but with additional oxygen in the tetra-II 
site, (d) as for (c) but with additional oxygon under the second 
Ag layer in the tetra-I and tetra-II sites. The average adsorp- 
tion energy with respect to the clean Ag(lll) substrate and 
free oxygen atoms, as well as the corresponding coverage, are 
given at the bottom of each figure. The relative variation of 
the first and second interlayer spacings with respect to the 
bulk value is also given to the right of the figures. Large pale 
grey and small dark circles represent silver and oxygen atoms, 
respectively. 
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at coverage 0.50, which as shown in Tab. I, is 2.57 eV for 
the latter and 2.94 eV for the former. 

We also tested if it is energetically more favorable for 
the two sub-surface O species to occupy different types 
of sub-surface sites. For instance, as shown in Tab. Ill, 

(Oocta+Otetra-Il)e=0.50 and (Ototra-I + Ototra-Il)6»=0.50, 

have average adsorption energies of 2.94 and 2.95 eV, re- 
spectively, very close to Oocta at same coverage (2.94 cV). 
We find that the energetically most favorable structures 
at coverage 0.75 ML do involve sub-surface O atoms in 
different sites, in addition to that where they both oc- 
cupy the octa site, namely, (Ohcp/Oocta + Otetra-ii)e=o.75 
with average adsorption energy 3.26 eV, (Ofcc/Otetra-i + 
Otctra-ii)e=o.75 with average binding energy 3.29 eV, as 
weh as, (Ohcp/Oocta)e=o.75 (3.25 eV). 

Interestingly, the slightly more favorable one, 
(Ofcc/Otetra-i + Otctra-ii)e=o.75 (depicted iu Fig. 10c), is 
similar in structure to the bulk silver oxide, Ag20, where 
both sub-surface oxygen atoms are located at the center 
of a silver tetrahedron, and the silver atoms shared by 
the O atoms are in a linear chain. Within the accuracy 
of our calculations, however, these three structures are 
essentially degenerate. 

B. Higher O concentrations and formation of an 
oxide-like film 

We now consider the incorporation of higher concen- 
trations of atomic oxygen. For a coverage of 0.75 ML the 
favored structures have 0.25 ML oxygen on the surface 
and 0.50 ML between the first and second Ag layers. We 
will now determine the energetic preference for adding 
one more O atom and then subsequently a second one, 
resulting in total coverages of 1.00 ML and 1.25 ML, re- 
spectively. In particular, we will investigate whether oxy- 
gen will prefer to occupy the sub-surface region under the 
first metal layer along with the two oxygen atoms already 
there, or if the O atoms will reside deeper in the "bulk" 
region, i.e., under the second metal layer. 

The results are shown in the lower part of Tab. III. For 
additional oxygen occupation under the first Ag layer, as 
an example, we show for the hcp/octa system the result 
for sub-surface coverages of 0.75 and 1.00 ML with O in 
the octa site. The average adsorption energies are 3.09 
and 2.89 eV, respectively. However, if these additional 
oxygen atoms occupy instead sites under the second Ag 
layer, in particular, according to the identified preference 
for a linear 0-Ag-O bonding arrangement, then the aver- 
age adsorption energies are 3.13 and 3.22 eV, respectively, 
which are clearly energetically more favorable compared 
to additional oxygen occupation under the first Ag layer. 
The latter structure is shown in Fig. lOd, where it can 
be seen, in comparison with 10b and 10c, that with con- 
tinuous O incorporation, the same structure can keep 
forming. Actually, the average binding energy of 3.22 
eV at coverage 1.25 ML, has already converged to the 



TABLE III. Average adsorption energy (per oxygen atom), 
iJad, and work funetion change, A$, for oxygen in various 

sites and for various coverages. 6'sub, ^on and ^totai represent 
the oxygen coverage in the sub-surface and on-surfaee regions, 
and the total coverage, respectively. The energy unit is eV. 





^on 


^sub 


Ea,d 


A3> 


^total = 0.25 










fee 


0.25 


- 


3.52 


1.23 


octa 


- 


0.25 


2.85 


0.09 


6'total =0.50 










fee 


0.50 


- 


2.92 


2.22 


hep 


0.50 


- 


2.83 


2.36 


fec-|-hcp 


0.50 




2.41 


2.53 


octa 


- 


0.50 


2.94 


0.09 


octa-|-tetra-II 


- 


0.50 


2.94 


0.02 


tetra-I-f-tetra-II 


- 


0.50 


2.95 


0.27 


octa-|-tetra-I 


■ 


0.50 


2.78 


0.46 


fcc/tetra-I 


0.25 


0.25 


3.27 


1.15 


hcp/octa 


0.25 


0.25 


3.33 


1.02 


6'total = 0.75 










nep/octa 


0.50 


0.25 


2.97 


2.22 


hcp/octa 


0.25 


0.50 


3.25 


1.42 


hep / octa-|-tetra-II 


0.25 


0.50 


3.26 


0.85 


icc/tetra-i 


0.50 


0.25 


2.90 


2.19 


iec/tetra-1 


0.25 


0.50 


3.06 


1.47 


fcc/tetra-I-|-tetra-II 


0.25 


0.50 


3.29 


0.98 


ftotal > 1.00 










hcp/octa 


0.25 


0.75 


3.09 


1.63 


hcp/octa 


0.25 


1.00 


2.89 


1.58 


fee /tetra-I-|-tetra-II 










/ tetra-I 


0.25 


0.75 


3.13 


0.66 


fee / tetra-I-)-tetra-II 










/ tetra-I-|-tetra-II 


0.25 


1.00 


3.22 


1.02 


fcc/tetra-I-|-tetra-II 










. . . / tetra-I-|-tetra-II 


0.25 


2.00 


3.21 


1.17 



value of the structure at the higher coverage of 2.25 ML 
(3.21 eV), as shown at the end of Tab. III. This value 
is also very close to the calculated heat of formation for 
bulk Ag20, which is 3.23 eV (also referred to the free O 
atom)."*^ 

In Fig. 11 the work function change of the identified 

energetically favorable structures involving on- and sub- 
surface oxygen, as well as that for the (4 x 4) structure are 
shown, where they are compared to the values for pure 
on-surface adsorption. It can be seen that they are sig- 
nificantly lower than for on-surfacc adsorption and vary 
between 1.1 to 0.5 eV. Interestingly these values arc close 
to that of --^LO eV as reported experimentally for Ag(lll) 
when subjected to high oxygen pressures and high tem- 
perature.^ The value for the (4 x 4) structure is also in 
very good agreement with the value of 0.55 eV reported 
by Bare et al. 

The identified energetically favorable structures with 
high O concentrations e.g. 1.25 ML or 2.25 ML are like 
silver (I) oxide Ag20(lll) but under a compressive lat- 
eral strain of 18%. The character of the oxide is already 
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FIG. 11. Work function change for on-surface O adsorption 
in the fee site (squares) (from Ref. 34) and of the (4 x 4) and 
the oxide-like structures with higher oxygen concentrations 
(triangles). 

present for coverage 0.75 ML: In addition to the same lo- 
cal atomic coordination, the vibration of sub-surface oxy- 
gen, Otcira-i, in (Ofcc/Otctra-i+Oictra-ii)9=o.7,>; IS Calcu- 
lated to be 526 cm~^ which is very similar to the exper- 
imental value of 545 cm~^ for oxygen in bulk Ag20.*^ 
The calculated valence band for the oxide-like struc- 
ture (Ofcc/Otctra-I + Otetra-Il/Otetra-I +Ototra-Il)e=1.25, 

shown in Fig. 12, also bears a close resemblance to that of 
bulk Ag20. The upper panel shows the result for the up- 
permost Ofcc-Ag-Otctra-i atoms, and the lower one shows 
the result for Ototra-ii under the first metal layer and 
Otetra-i Under the second metal layer, as well as the sil- 
ver atom in the second layer to which they commonly 
bond. Compared to the clean surface, or lower coverage 
structures, e.g., (Ofcc/Otetra-i)e=o.50 as shown in Fig. 7, 
the silver 4rf-band has considerably narrowed due to the 
significant increase of the spacing between the metal lay- 
ers. Also, two extended shoulders appear either side of it. 
These features are very similar to those; of the the oxide 
surface under zero lateral strain, which will be reported in 
detail elsewhere.^® We calculated the removal energies of 
the on-surfacc oxygen of the favorable structures listed 
in Tab. III. The resulting energies are very similar to 
what we found for the favorable 0.50 ML on-surface-|-sub- 
surfacc structiircs discussed above in Sec. IV, namely, 
they are all in the range of 3.90~4.01 eV. These bond- 
strengths are actually the strongest we have found for 
this system and such species could be important for un- 
derstanding how silver functions as an efficient oxidation 
catalyst. 

1. Relative stability of the oxide-like structures 

It has been seen above that the average O adsorp- 
tion energies of the energetically favorable structures at 
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FIG. 12. Total projected density of states for the oxide-like 

structure 

(Ofcc/Otctrai+Otctra-Il/Otetrai-|-Otetra-n)e = 1.25. The UppCr 

figure shows the PDOS for Ofcc (dashed line), the uppermost 
Ag atoms bonded to both O atoms (solid line) and Otetra-i 
(dotted line) under the first layer. The lower figure shows the 
PDOS for Otetra-ii (dashed line) under the first metal layer, 
the shared Ag atom (solid line) and Otetra-i (dotted line) un- 
der the second metal layer. The Fermi energy is indicated by 
the vertical dashed line. 

0=0.50 and 0.75 ML are slightly greater than the heat 

of formation of the bulk oxide referred to atomic oxy- 
gen (compare 3.27 and 3.29 eV to 3.23 eV). For the 
higher concentration structures the values are however 
extremely close (compare 3.22 and 3.21 to 3.23 eV). 
This is very interesting since the oxide- like structures, al- 
though similar in geometry to ideal Ag20(lll) surfaces, 
are laterally compressed by a large amount. To inves- 
tigate why the energies are so close, e.g., whether the 
underlying Ag(lll) surface plays a role, or if the distor- 
tion actually has little effect on the energetics, we study 
the relative stability of the oxide-like structures in com- 
parison to true Ag20(lll) structures. In particular, wc 
calculate free-standing Ag20(lll) layer structures un- 
der zero strain and under a compressive strain of 18%, 
which corresponds to the strain of our identified oxide- 
like structures due to being commensurate with the un- 
derlying Ag(lll) surface. 

For a thickness involving 5 metal layers, it is found that 
the true oxide is energetically less favorable by 0.05 eV 
per oxygen atom than the oxide film under the strain. 
For thicker films corresponding to 7 metal layers, the en- 
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crgy difference between two structures decreases to just 
0.005 eV per oxygen atom, where the deformed oxide film 
is still preferred. Only for thicker structures involving 9 
metal layers does the true oxide film (zero strain) become 
energetically favorable (by 0.02 eV per oxygen atom). 
This could be taken to indicate that a transition to the 
natural oxide surface will only happen at or beyond a 
critical thickness corresponding to about nine metal lay- 
ers. The transition may not occur immediately because 
the energy gain is still very small and it may take more 
layers in order for sufficient strain energy to build-up. 
Also, the interaction between the oxide structures and 
the metal substrate, which was not considered, may hin- 
der the transition. That the structures with the smaller 
lateral dimensions are not notably less favorable than the 
true oxide structures can be understood in that there is 
a significant expansion normal to the surface, as shown 
in Fig. 10. These energetics are in contrast to what was 
found for the O/Ru(0001) system,^" where from similar 
calculations it was foimd that RuO2(110) layer structures 
become energetically more favorable compared to layers 
of the identified on-surface-|-sub-surface structures with 
the same stoichiometry, after a thickness equivalent to 
only two or more metal layers. In this case the lateral 
deviation between the natural RuO2(110) oxide and the 
on-surface-l-sub-surface structure was equal to 15 % and 
35 % in the y and x directions, respectively, where the 
natural oxide is more expanded. 

2. Adsorption of an ozone-like species 
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FIG. 13. Total projected oxygen density of states for tfie 
ozone-like species adsorbed at an Ag surface vacancy. The 
upper panel shows the valence band energy region and the 
bottom panel, the lower energy region. "Om" indicates the 
center O atom of the ozone- like species and the "Ol/r" indi- 
cate the left and right O atoms of the ozone-like species that 
are bonded to two surface Ag atoms. 



Another issue of debate in the literature is the ori- 
gin of several electronic states observed lying in the en- 
ergy region below the Kg-Ad band and above the 0-2s 
semi-core level. On one hand it has been proposed that 
hydroxl groups are responsible^ and on the other, an as- 
sociatively adsorbed ozone-like species at a surface Ag 
vacancy has been proposed. We investigated the latter 
species, the atomic geometry of which is shown in the 
insets of Figs. 13 and 14. We find, in agreement with 
Ref. 50, that the molecular ozone-like species is ener- 
getically preferred compared to three separate O atoms 
adsorbed at the vacancy edges. We considered two va- 
cancy concentrations for the adsorption of the ozone-like 
species, namely, 0.11 and 0.33 ML, corresponding to O 
coverages of 0.33 and 1.0 ML, respectively. The atomic 
geometry is found to be similar to the ozone molecule 
in the gas phase, namely, when adsorbed on the surface 
the calculated bondlength between the O atoms is 1.41 A 
and the angle is 118 degrees, while for the free molecide 
the respective values are 1.27 A and 117 degrees. The 
longer bondlength between the O atoms when adsorbed 
on the surface indicates that the inter 0-0 bonding is 
weakened. The adsorption energy is very similar in each 
case, indicating that they interact little with each other. 
Compared to the oxide-like structures they are energet- 



ically unfavorable. However, if vacancies were provided 
'for free' (i.e. not taking into account the vacancy forma- 
tion energy), then such a species is energetically favorable 
for 6 = 1.0 ML (see Fig. 14). We expect, however, that 
the concentrations of vacancies will be rather low and 
this species will not be present in large concentrations. 

With regard to the electronic structure, we find that 
there are indeed a number of electronic states in the en- 
ergy regime below the bottom Ag-4d band as can be seen 
from Fig. 13, which shows the projected density of states 
on the O atoms. In particular, features are seen at ener- 
gies of ~2.0, 3.5, 6.0, 7.0, 8.5, 14.4, 19.5, 24.8 eV below 
Ep. Reported experimental features observed on poly- 
crystalline silver are at 3.5, 9.2, and 11.2 eV below Ep, 
which were thought to be due to hydroxls,^ and features 
observed for the O/Ag(110) system are at 3.3, 7.0, 9.5, 
12.7, 17.3 eV, which were attribued to the associatively 
adsorbed ozone-like species. It can be appreciated that 
there is some overlap of these various energies but it is 
not possible to draw any solid conclusions on the basis of 
such comparisons. 

In order to visually summarize the main energetics 
of our results, we show in Fig. 14 the average adsorp- 
tion of oxygen at Ag(lll) versus coverage. This also in- 
cludes our earlier results for pure on-surface and surface- 
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substitutional adsorption^^ as well as for investigations 
into bulk substitutional adsorption. There is one point 
in this figure not discussed in the present or previous 
(Ref. 34) paper so far and that is for the adsorption of 
an oxygen atom on the surface in the fee site near a pre- 
existing vacancy with coverage 0.25 ML. The adsorption 
energy is seen to be less favorable than on the ideal (111) 
terrace, but if such vacancies were provided "for free" 
(i.e. without paying the vacancy formation energy cost) 
it can be seen that their binding to the surface is con- 
siderably stronger than on the ideal (111) terrace, com- 
pare 3.83 eV to 3.52 eV. We expect a similarly stronger 
binding of O at other under-coordinated metal sites such 
as at the edge of steps or kinks, as has been found for 
other 0/transition- metal systems. Considering Fig. 14, 
at a glance, these results predict that at low coverage, 
O prefers to adsorb on the surface, but with increasing 
coverage the thin (4 x 4) surface-oxide is favored. For 
higher concentrations, thicker oxide-like structures simi- 
lar to the (111) surface of Ag2 are preferred which have 
a coverage of 0.25 ML on the surface and 0.50 ML be- 
tween each Ag layer. These energetics relate well to the 
general behavior found experimentally in the early work 
of Czanderna,^^ whose study shows that the isosteric heat 
of adsorption of oxygen on silver powder decreases from 
0.91 to 0.37 eV per oxygen atom (with respect to 1/2 
the binding energy of O2), and then remains constant at 
0.40±0.02 eV per oxygen atom for coverages 0=0.33 to 
about 0.90 ML, which was attributed to the growth of an 
oxide layer. 



C. Thermodynamic and kinetic considerations 

For a metal in contact with an oxygen atmosphere, in 
thermal equilibrium, the bulk oxide will form when the 
oxygen chemical potential, /io, satisfies the equation. 



MO + 2/XAg(bulk) — /^AgaO 



(4) 



where /XAg(buik) and /iAg2 are the chemical potentials 
of an Ag atom in bulk and of the bulk oxide, Ag20, 
respectively. This equation can be expressed in terms of 
the total energies of the respective systems, i.e.. 



Er 



2E 



Ag 



■^Ag2 



(5) 



Using the above equation, and referring the energy to 
that of a free O atom, we obtain for the O chemical po- 
tential Ho = 3.23 eV. This is just the heat of formation 
of Ag20, AiJoxide, but referred to a free O atom instead 
of to I/2O2. 

Concerning the transition from chemisorption to bulk 
oxide formation, it has been proposed that the critical 
coverage, Oc, at which the transition occurs is thermody- 
namically, rather than kinetically, determined. Below 9c, 
the differential heat of adsorption of chemisorbed oxy- 
gen is higher than the heat of formation of the oxide. 
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FIG. 14. Average adsorption energy (with respect to 
atomic O) versus coverage for the different structures con- 
sidered: pure on-surface oxygen in fcc-hollow sites (dia- 
monds), surface- and bulk-substitutional adsorption (trian- 
gles), pure sub-surface oxygen in octahedral sites under the 
first (subsurface-octa.) and second (bulk-octa.) Ag layers 
(circles), an ozone-like molecule adsorbed at a surface Ag va- 
cancy (left-pointing triangles), and surface oxide-like struc- 
tures (squares), as illustrated in Fig. 10. The vertical down- 
ward pointing arrows directed at dashed lines with open sym- 
bols, indicate the corresponding binding energies of the struc- 
tures involving surface Ag vacancies, i.e. the vacancy forma- 
tion energy, i?™'^ (see Ref. 34 for the definition), is not taken 
into account. The adsorption energy of O in the fee site (tri- 
angle) with a neighboring surface Ag vacancy is also shown 
for 9 — 0.25 ML, as well as the corresponding binding energy. 



but above 6c it becomes lower. The differential heat of 
adsorption at a given coverage, H^^{9), can be thought 
of as the energy required to remove a single O atom on 
the surface at that coverage, or equivalently, the energy 
gain on adsorption of a single O atom to result in that 
coverage (i.e., with the presence of other pre-adsorbed 
O atoms). The criteria for the critical coverage beyond 
which bulk oxide formation is expected to occur is. 
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(6) 
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Applying this to the 0/Ag(lll) system, using the on- 
surface adsorption energies from our previous paper, '^^ 
we obtain a critical coverage 6c of 0.24 ML. This value 
correlates very closely to the coverage of 0.25 ML found 
from our extensive DFT calculations when oxygen starts 
to penetrate into the sub-surface region, suggesting onset 
of oxide formation and occupation of sub-surface sites 
are linked. Similar correlations for O at other transition 
metal surfaces, namely, Ru(OOOl), Rh(lll), and Pd(lll) 
occur, where it is found that the critical coverage is higher 
for the "harder" metals towards the left of the TM series 
(i.e., Ru).23 

In contrast to our theoretical results which predict the 
formation of oxide-like films, to our knowledge, exper- 
imental studies employing high pressures and elevated 
temperatures have only reported the thin (4 x 4) surface- 
oxide to date.''^''''''^^'^^'^^ Formation of an oxide requires 
firstly sufficient numbers of dissociated oxygen atoms, 
and secondly, it requires significant mass transport (e.g., 
diffusion of Ag atoms or O atoms) which can only take 
place at elevated temperatures. Given that the heat of 
formation of bulk silver oxide is low, we expect that for- 
mation of thicker oxide-like layers may be prevented ei- 
ther because the temperature at which they likely decom- 
pose is not high enough for sufficient atomic rearrange- 
ments to occur, or that so far there have been insuf- 
ficient concentrations of oxygen atoms delivered to the 
surface, indeed the (111) surface of silver has an ex- 
tremely low dissociative sticking coefficient for O2, 10~^. 
Possibly with the use of atomic oxygen or ozone (which 
readily dissociates at the surface) and low temperatures, 
thicker oxide films could be formed experimentally. This 
has recently been done for studies on the oxidation of 
the Ft (111) surface^^ as well as for polycrystallinc sil- 
ver when exposed at room temperature to a mixture of 
oxygen and ozone, where thick films have recently been 
reported.^** 

In order to gain more insight into the behavior of the 
O / Ag system as an oxidation catalyst and to identify pos- 
sibly active O species, it is important to be able to take 
into account the effect of temperature and gas pressure. 
In particular, the catalytic reactions of ethylene cpoxida- 
tion and the partial oxidation of methanol take place at 
two significantly different temperature ranges (500-600 K 
and 700-900 K, respectively) and both occur at high (at- 
mospheric) pressure. Our investigations in this direction 
are in progress. 

VI. CONCLUSIONS 

In the present work, we have systematically investi- 
gated pure sub-surface oxygen and structures involving 
both oxygen adsorbed on the surface and sub-surface 
oxygen, as well as oxide- like structures at the Ag(lll) 
surface through density functional theory calculations. 
Compared to pure on-surface oxygen adsorption, which 



exhibits a strong decrease in adsorption energy with in- 
creasing coverage, pure sub-surface adsorption exhibits a 
weak dependence on the coverage, where the octahedral 
site is favored for all of the investigated coverage range of 
0.11 to 1 ML. Comparing the adsorption energies of pure 
on-surface O and pure sub-surface O, we find the former 
is energetically favored over the latter for coverages up to 
around 1/2 a monolayer, whereafter the pure sub-surface 
structures arc preferred. 

The energetically favorable structures involving both 
on-surface and sub-surface oxygen have a lower density 
of states at the Fermi energy and involve less ionic (less 
negatively charged) O atoms as compared to the energet- 
ically unfavorable structures. As quantified by the con- 
cept of a "removal" energy, the presence of sub-surface 
oxygen modifies the on-surface oxygen-silver bond sig- 
nificantly. Depending on the adsorption site, it can ei- 
ther stabilize or destabilize on-surface oxygen, and vice 
versa. The energetically most favorable structures, how- 
ever, stabilize on-surface oxygen. 

On the basis of the energetics of all the calculated 
structures we find the following scenario: At low cov- 
erages oxygen prefers to stay on the surface in the fee 
site. With increasing coverage a thin (4 x 4) surface-oxide 
is energetically favorable, and for higher concentrations, 
the calculations predict the formation of oxide-like struc- 
tures very similar to Ag20(lll). We propose that with 
the use of atomic oxygen or ozone and low temperatures, 
such thicker oxide-like films may be observed experimen- 
tally. 

Having thoroughly studied the energetics, and atomic 
and electronic properties of the 0/Ag(lll) system, as is 
crucial for gaining further understanding into the func- 
tion of silver as an oxidation catalyst, the next consider- 
ations must address the effect of elevated temperatures 
and high pressures, since these are the conditions under 
which real catalysis takes place. This will be the subject 
of our following publications. 
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